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Abstract 

We estimate the flux of garnma-rays that result, from collisions of high energy galac- 
tic cosmic rays with the solar atmosphere. An important aspect of our model is the 
propagation of cosmic rays through the magnetic fields of the inner solar system. We 
use diffusion to model propagation down to the bottom of the corona. Below the corona 
we trace particle orbits through the photospheric fields to determine the location of cos- 
mic ray interactions in the solar atmosphere and evolve the resultant cascades. For our 
nominal choice of parameters, we predict an integrated flux of gamma rays (at 1 AU) of 
& (E 7 ^ 100 Me\ ) ^ 5 x 10 ^ cm " sec F This can be an order of magnitude above the 
galactic background, and should be observable by EGRET. 


Int roduction 

There aie a numbei of reasons to study the sun as a source of 7* rays produced in cosmic rav 
cascades. First, there is the purely intellectual question of whether the 7-ray albedo is strong 
enough to be observed. Second, there is the practical issue of whether the sun could be a 
strong enough source of 7 s to be a confusing background to other sources. Third, many of 
the same issues that occur in studying solar flares as 7-ray sources are likely to come into the 
calculation of the 7-ray albedo from cosmic ray interactions. These include the structure of the 
magnetic fields near the sun and charged particle propagation through those fields. The details 
of those issues will be different in the two cases, and so studying the 7-ray albedo may give one 
a new perspective on solar flare models. Finally, whereas models of other astrophvsical 7-ray 
sources may be largely speculation, due to the sun’s proximity models of the solar source may 
be directly tested in other ways. Successful modeling of the sun may lead to better models of 
other more remote sources. 

In this contribution, we summarize a recent paper 1 in which we discuss a model for esti- 
mating the albedo 7-ray flux from cosmic ray cascades in the solar atmosphere. Our model 
explicitly includes several Tieliomagnetic’ effects. A naive estimate might assume that the cos- 
mic iay flux incident on the sun s surface equals the flux at Earth, and that the solar albedo, 
like the Earth’s, might be quite small. We argue that the situation is just the reverse: inter- 
planetary (IMF) and coronal magnetic fields suppress the flux reaching the solar surface, but 
photospheric magnetic fields result in an efficient albedo. These conclusions coupled with a 
Monte Carlo study of the photon production by cascades in the solar atmosphere lead to a 
prediction of a flux that is detectable by the EGRET 10 instrument on the Arthur H. Compton 
Gamma Ray Observatory (GRO). We raise and disc uss some questions for future study. 

Cosmic Ray Absorption by the Sun 

framework for Absorption. We write the absorption rate for cosmic rays of energy E as 

V(E) = 47 rr 2 j(E), (1) 
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whore } is the net radial (lux of cosmic rays. For a simple model ignoring magnetic fields. 

charged particles travel in straight lines and j is given b> j/v ( -V ) A . 

tTf) is the differential density of cosmic rays alra oo, ,9c is their velocity and the subscnp 
N indicates that this is the naive value for absorption on a sphere of radius *15. We 
/<( = M E), the observed density at Earth 2 . The naive absorption rate is then 


r,v - 


More generally, we write the absorption rate as 

ACq\'n 


1 ' 


A + Cp - AC 


( 2 ) 


(3) 


D 


where Co is a correction factor to account for propagation through the IMF and corona, and A 
is the probability that a primary which reaches the photosphere will be absorbed there ms e 
"f heing“d back into the corona. The denominator arises from considering the possibility 
that a cosmic ray which was reflected back into the corona may diffuse back tot lie pho osp 
and be absorbed on the 2" (< , 3 rd , . . . penetration of the photosphere before totally escap ng 

solar system. 

Diffusion. We estimate C D using a diffusion model for cosmic ray Propagation in thejne, - solar 
system Since the IMF is expected to be nearly radial inside the Earth s orbit *e take 
diffusion problem to be spherically symmetric with a partially absorbing boundary condition 

alte pliotosphere. The inward flux is given by j = - Df, my be used ,n h,. 1 to solve 

for r. For the diffusion ‘constant’ we consider the form, 


D(E, r) = % (r/r e )° (E/\ GeVf 


(4) 


where D is normalized to D„ its value at K = 1 CieV and r = r 4 Eg. 1 is subject to the 
constraint DIE rl < rc/8, which is dictated by causality, i.c. the net flux j cannot, excc JN- 
Z nominal choi e of diffusion parameters 2 is = (0.09. 2. 1), where ,s given in 

uni s o " X A u. Unfortunately, the diffusion parameters are mostly based upon measurements 
bv'spacerraif in Earth orbit or the outer solar system, and this leads to uncertainty in our 
predictions. On the other hand, an observation of the signals we propose may provi e a ne 
way to probe the inner solar system. 

The Absorption Coefficient, A. To determine the absorption probability we take an isotropic 
flux of cosmic rays at the bottom of the corona and propagate them inward including 
magnetic effects and absorption by the ambient gas. For this part of our calculation we aban 
diffusion and trace the particle orbits exactly. The small scale magnetic structure - used for th 
calculation is illustrated in Fig. 1. It consists of a loose network of flux tubes which penetra.. 
the photosphere at the corners of convective cells, where the field is swept by the moving u . 
The tubes open up to a space filling ‘canopy’ at the bottom of the corona wh.ch we take. to 
be 1250 km above the photosphere. We use field strengths within the tubes of 1500 G at 
the ^photosphere and 6.5 G at the bottom of the corona. We assume the flux tubes maintain 
pressure equilibrium with the gas, t.e. B 2 ~ P„. Due to the pinching of the flux Uibes, most 

primaries are mirrored and return to the corona without interacting, iowc ' 

oss cone where nearly vertical cosmic rays penetrate to sufficient depth to be absorbed 
the nominal parameters we find an absorption probability of A = -0052 for protons; and .0088 

for 4 He. 
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Figure 1: Model of magnetic fields near the photosphere. Shading increases with 
magnetic field intensity. 

The discussion so far is valid only for relatively low energy primaries, which we argue will be 
trapped on field lines and necessarily enter the flux tubes. Higher energy primaries will not be 
trapped if their gyroradii are comparable to the distance bewteen flux tubes. We estimate the 
threshold for trapping to be Ej ^ 3 TeV, and take A(E > Ej) = 1. As cosmic rays with such 
high energies are relatively rare the uncertainty in Ej should not strongly affect our results. 

Absorption Rates. We show our nominal absorption rate for protons as the bold curve in 
Fig. 2. We also show, as dashed lines, the results for a pessimistic (D e ,a,/3) = (0.01, 2, 0.5), 
and optimistic (0.1, 1.5, 2) choice of diffusion parameters; and, as light, solid curves, three 
different absorption coefficients (A ~ 0.001, 0.01, 0.1). We also show the naive flux from Eq. 2. 
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Figure 2: Absorption rate of cosmic ray protons by the sun. The bold curve is 
our nominal result. Varying the diffusion parameters results in the dashed curves. 
Varying the absorption coefficient yields the light, solid lines. The naive result 
(Eq. 2) is shown as a dotted curve for comparison. 
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The 7-ray Albedo 

In calculating the 7 -ray albedo we distinguish between two energy regimes; E < E T and 
E > E t . As stated above, for E > E T , cosmic rays are not trapped on flux lines and we 
may ignore local magnetic flelds. The 7 -ray albedo of these cascades will be very small, with 
contributions only from incident protons grazing the photosphere. Primaries with E < T, 
on the other hand, are assumed to stay within the flux tubes. Due to magnetic mirroring t e 
albedo from these cascades may be significant. 

Two Models for Cascade Propagation. Nearly vertical cosmic rays in the interior of the loss 
cone interact before being mirrored, Le. while still having a downward velocity ( lg. a). e 

boundary of the loss cone, however, is a region where protons may interact after being mirrore 
(Fig. 3b). These interactions will contribute to an upward flux of gamma-rays. We find that 
there is roughly a 25% chance that primary interactions will take place after the primary has 
been mirrored. Since 7 -ray mean free paths are comparable to those for the primary hadrons, 
the albedo of the sun is considerable. 



Figure 3: Cascade geometries for production of the 7 -ray albedo. Shaded circles 
indicate interaction sites. Heavy lines show the primary trajectory before the 
interaction. The curved arrow in a) indicates the path the primary would take it 
no interaction occurred, a) Interaction occurs before the primary is mirrored. No 
albedo is produced if the cascade develops linearly, b) Interaction occurs after the 
primary is mirrored. This case contributes to the albedo, c) Interaction occurs 
early, but if the whole cascade is mirrored albedo photons may still result. 

Another aspect of the low energy cascades is that the charged particles will follow the same 
trajectory the parent cosmic ray would have followed if it had not interacted; r.e. the whole 
cascade is to some extent mirrored (Fig. 3c), even though some of the cascade energy is neutral. 
To explore this effect we took two extreme models; one where the cascade developed along a 
constant zenith angle, and one where it developed along the path of a charged particle being 
mirrored out of the photosphere. In the first model photons can arise only from incident cosmic 
rays that interact after being mirrored, but in the second even cosmic rays which interact while 

moving downward can produce an upward flux of 7 -iays. 

Given these preliminaries, we calculate the 7 -ray albedo by evolving ‘atmospheric cascades 
from the point where the primary cosmic ray interacts, through the solar material, to the 
surface. In the first model, the path length distribution is given by the slant depth from the 
interaction site. Cascades from primaries which interact before being mirrored are given infinite 
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path lengths and produce no albedo 7-rays. In the second model the path lengths are computed 
as if all particles are charged, in which case all cascades may in principle produce albedo 7-rays. 

7-ray Production. The 7-ray flux at Ka.rtli is given by 




1 dLy 
1 7T 7'^, (I " 


where the 7-ray luminosity of the sun is given by 



(• r >) 


( 6 ) 


The sum over i includes both p and 'lie primaries, and l\ is given by Kq. 3. 1 he 7-ray yield 
per absorbed primary of energy E is given by averaging over the path length distribution 

r d P 

V, ( E y , Ei ) = J y t ( E y . E , , . 7 - ) <1 x , ( 7 ) 

where y, is the yield of photons of energy E~ t that result, from cascades with path length x 
initiated by primaries of energy E, , and dP t /dx is the probability that a primary of species i 
produces a cascade with path length x. In separating /’ ( .r ) and 1’ , ( /■' ) we have ignored the slow 
energy dependence of the hadronic cross-section. The yields, y,. are calculated using standard 
Monte Carlo techniques 6,1 including both hadronic and electromagnetic interactions. Note that 
the density is low enough that for primaries with E < Ej, all mesons decay before they are 
absorbed or stop. 



Log E y (MeV) 


Figure 4: Photon flux at lOarth. The solid (dashed) curves show results for the 
nominal (naive) absorption rate. In each pair, the upper curve shows a model 
where the cascades are mirrored, while the lower curve assumes the cascades 
follow straight trajectories. The heavy curve shows the galactic background at 
b = 30°. 
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Results The results of our 7 -ray calculations are illustrated in Fig. 4. The expected background 
is due to cosmic ray collisions with intergalactic gas 8 , combined with observational data be ow 
200 MeV, ami i, S iiow„ for galactic, latitude I, = 110". The background shou Id scale *6 1 /sm ». 
Boll, signal and background arc given for a solid angle equal to the solar disk. The . g 
larger than the galactic background for reasonable assumptions about cosmic ray absorption. 

Other Albedo Products 

We have performed similar calculations 1 to estimate the albedos for other cascade products, 
which we summarize here. 

Neutrons. From the perspective of the GRO the most interesting possibility is a quiet time 
flux of high energy neutrons from the sun. The primary production mechanisms for neutrons 
are spallation reactions of 4 He and charge exchange interactions between me, dent proto: is and 
target hadrons. Neutrons are unstable, so there is no galactic background to speak of but as 
the distance to the sun is comparable to the neutron lifetime, neutrons produced on the 
may still reach the Earth before they decay. 

The geometry of the neutron production mechanism is similar to that for photons - the 
albedo arises from cascades that are mirrored or for which the primary was mirrored before 
interacting. For neutron production we modified our Monte Carlo of the cascades to include 

explicit charge separation of protons and neutrons. 

We estimate a quiet time neutron flux of I n (E n > 100 MeV) .5 X < 1 * cm ^ 

approximately a factor of 10 4 below the pre-GRO observational limits . Both COMP ' 
and OSSE 14 have sensitivity to neutrons, as evidenced by their detection of neutrons during 
the June 1991 flares. However, their capabilities are not nearly good enough to detect our 
predicted quiet time flux. For example, COMPTEL’s source sensitivity to photons is not quite 
1(T 4 cur 2 sec" 1 , and it’s neutron sensitivity is even less. It is interesting to note tha ie 
integrated neutron flux during the June 9 flare corresponds to ~ 30 yr of our predicted qme 
time flux, albeit in a somewhat lower energy band. 

Neutrinos. Neutrinos produced by solar cascades are potentially observable through their 
conversion to a flux of upward going underground muons. The calculation for neu nnos i er 
qualitatively from that for photons in that the high energy cascades dominate the s.gna i Th 
relative importance of the high energy cascades arises since they are essentially unaffected by 
diffusion in the IMF and corona and because there is a roughly lineal increase in >o i ie m 
production cross-section and range with E v . Also note that neutrinos (unlike photons) are not 

generally absorbed bv the sun, and so there is no cutoff at h > Er- 

For the neutrino signal the principle background is the flux of neutrinos produced in ter- 
restrial cascades 11 . For cascades with E > 500 GeV the sun is ~ 10 fames more efficient at 
producing neutrinos than the Earth, due primarily to the larger sea e height of the solar at- 
mosphere. As a result if a neutrino telescope could resolve the sun, detection would be sign 
limited - not background limited. However, for our nominal model we calculate an undergroun 
muon flux of /( E„ > 10 GeV) ~ 5 X 10 - 17 cm ' 2 sec" 1 , which is less than one event per year 
at any operational neutrino telescope, and so we conclude that the neutrino signal is curren y 

unobservable. 

Antimatter. Finally, we have considered antimatter signals, but do not find them promising. 
Antiprotons and positrons are charged, and so lack a directional signal. Antineutrons would be 


547 


difficult to distinguish from a much more copious neutron flux. 

Future Directions 

We have made a reasonable estimate of the 7 -ray albedo from cosmic ray absorption by the 
sun; however, there are a number of issues that still need to be addressed. 

Observing Strategies and Phase I EGRET Observations. The first question - is the signal 
really there? If the sun were resolved by EGRET then the source strength would be many 
times the galactic background. However, EGRET’s point spread function 6 is greater than the 
size of the sun 15 (b = 1.4° for 150 MeV < E < 500 MeV) and so it is important to reduce 
the background by observing the sun while it is at high galactic latitude. During the phase 1 
observing program there are 29 scheduled 2 week observing periods (and several shorter ones), 
of which 4 (fortuitously) have the sun within EGRET’s 20 ° 11WHM field of view during part 
or all of the observing period. The three best of these (Oct. 3, 1991; Aug. 6 , 1992; and April 
2, 1992) have the sun’s galactic latitude at | 6 | > 44°. It is therefore reasonable to estimate 
EG RE 1 s exposure during the phase one viewing period as ~ 1.8 X 10 9 cm 2 sec, which would 
give a yield of approximately 100 photons with E > 100 MeV. By comparison, the total 
number of galactic background photons detected within S of the sun during the same exposure 
is expected to be — 240. A three sigma excess would then correspond to roughly 50 photons, 
and assuming our estimates of the signal are correct, the sun should be an observable source. 

Improvements to the Calculations. In the event of a successful detection, there are a number 
of improvements that could be made without abandoning the basic structure of our model. 
We anticipate that some of these would increase our predicted flux if included and some would 
decrease it, but taken as a. whole they would not change the results bv an amount large compared 
to the uncertainty in a particular choice of diffusion and canopy parameters. 

In no particular order - a) Eliminate the discontinuous treatment of absorption at Ej\ This 
should not dramatically affect the flux estimates for E 1 in the range 100 MeV - 1 GeV since 
their parent primaries typically have E «C Eq\ b) Include an energy dependent leakage term in 
the calculat ion of the absorption rate. In the present model, a low energy primary which reaches 
the sun is either absorbed in a flux tube or eventually diffuses back out of the solar system. 
It is also possible that some fraction will leak across field lines to be absorbed in the space 
between flux tubes. Since these primaries and the resultant cascades would not be mirrored 
their contribution to the albedo would be suppressed, c) In the present Monte Carlo, charged 
particle secondaries that escape the sun were thrown away; however, it would be consistent with 
the treatment of non-absorbed primaries to assume that these charged secondaries eventually 
return to the sun and are absorbed, continuing the cascade after a ‘pause 1 . We estimate that 
this would increase the 7 -ray signal by a factor of 2-3. d) Include charge separation in the 
electromagnetic part of the cascade, e) Make the model of the flux tubes two dimensional 
to account for their finite radial extent, and to allow the magnetic field to have a horizontal 
component. J) Incorporate modulation by the solar wind into the incident flux calculations. 

Phenomenological Considerations. Assuming that the signal is seen, there are various features of 
the 7 -ray signal that may be useful either as diagnostics of the absorption model, or as indicators 
of transport phenomena around the sun. An obvious example is measuring the spectrum of the 
7 -ray flux and comparing it to the model predictions. Another possibility arises from noting 
that most of the interplanetary magnetic field is anchored to one of the polar regions of the 
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sun If the primaries find it difficult to cross field lines (i.e. diffusion takes place parallel to 
bnt not perpendicular to the local magnetic held) then one might expect the polar caps to be 
ntaiffclntTy brighter in p-mys than the etpiator. Although the diameter of the sun ,s smaller 
than EGRET’S point spread function one might still hope to pick out a wea . tpo e i pa , 
upon gathering sufficient data, especially at higher photon energ.es where EGRF/1 

" ^Th'is not a complete list, but as these two examples indicate a study of the 7 - ray flux 
from the sun can potentially provide useful information about 7-ray pro , uction mec a" is 
and cosmic ray transport in the solar system, and by extension to other sites m 
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